Non-Boussinesq Convection at Low Prandtl Numbers: Hexagons and Spiral Defect 

Chaos 
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We study the stabihty and dynamics of non-Boussinesq convection in pure gases (CO2 and SFg) 
with Prandtl numbers near Pr ~ 1 and in a H2-Xe mixture with Pr — 0.17. Focusing on the 
strongly nonlinear regime we employ Galerkin stability analyses and direct numerical simulations of 
the Navier-Stokes equations. For Pr ~ 1 and intermediate non-Boussinesq effects we find reentrance 
of stable hexagons as the Rayleigh number is increased. For stronger non-Boussinesq effects the 
hexagons do not exhibit any amplitude instability to rolls. Seemingly, this result contradicts the 
experimentally observed transition from hexagons to rolls. We resolve this discrepancy by including 
the effect of the lateral walls. Non-Boussinesq effects modify the spiral defect chaos observed for 
larger Rayleigh numbers. For convection in SFe we find that non-Boussinesq effects strongly increase 
the number of small, compact convection cells and with it enhance the cellular character of the 
patterns. In H2-Xe, closer to threshold, we find instead an enhanced tendency toward roll-like 
structures. In both cases the number of spirals and of target-like components is reduced. We quantify 
these effects using recently developed diagnostics of the geometric properties of the patterns. 

PACS numbers; 47.20.Bp, 47.52. +j,5.45.Jn,47.54.-|-r 



I. INTRODUCTION 

Rayleigh-Benard convection of a horizontal fluid layer 
heated from below is a paradigmatic system for the study 
of complex spatial and spatio-temporal patterns Q . Usu- 
ally, the system parameters are chosen such that vari- 
ations of the fluid properties across the layer are very 
small. In the theoretical description this allows the 
Oberbeck-Boussinesq (OB) approximation, in which the 
temperature dependence is neglected in all fluid proper- 
ties except for the density term that is responsible for 
the buoyancy. For large temperature differences across 
the layer the variation of the fluid properties with the 
temperature become significant. These non-Oberbeck- 
Boussinesq (NB) effects break the up-down symmetry 
that is characteristic of the OB approximation and, thus, 
allow a resonant triad interaction among three Fourier 
modes whose wavevectors form a hexagonal planform. 
Due to the resonant interaction the primary bifurcation 
to the hexagons is transcritical, and hexagons arepre- 
ferred over rolls in the immediate vicinity of onset |3l . 

Within the framework of the leading-order weakly 
nonlinear analysis, hexagons typically become unstable 
to rolls further above threshold, where the amplitudes 
are larger and the resonant-triad interaction loses sig- 
nificance coinpared to interactions involving four modes 
0, 0, 0, 0, 0, 0- This scenario of a transition from 
hexagons to rolls has been confirmed in c[uite a number 
of experimental investigations 1^, ^, 10^ 0, ^ , and 
quite commonly it has been assumed that hexagon pat- 
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terns in NB convection are confined to the regime close 
to onset. 

Two convection experiments using SFg as the work- 
ing fluid MM 

have shown, however, that even in the 
strongly nonlinear regime stable hexagon patterns can be 
observed. Under OB conditions hexagons were found 
at relatively high Rayleigh numbers, e = [R — Rc)/Rc~ 
3.5 Due to the mid-plane symmetry of OB convec- 
tion hexagons with up-flow in the center coexist stably 
with down-flow hexagons in this regime. In experiments 
using SFg under NB conditions [13 the hexagons that 
appear at threshold were replaced by rolls for somewhat 
larger heating and then reappeared in the strongly non- 
linear regime near e = 0(1). The restabilization was 
attributed to the large compressibility of SFg near its 
critical point T^. The hexagons that regain stability 
were termed reentrant hexagons. 

Recent numerical computations 16] have demon- 
strated that hexagons can restabilize in NB convection 
even if the fluid is incompressible. For instance, in water 
hexagons that are unstable at e = (-R — Rc)/Rc = 0.15 
can already restabilize at e = 0.2. The origin of this reen- 
trance was traced back to the existence of stable hexagons 
in OB convection at large Rayleigh numbers, and the de- 
pendence of the NB effects on the Rayleigh number. 

At low Prandtl numbers (Pr ~ 1) and further above 
onset OB convection exhibits a new state: spiral de- 
fect chaos (SDC). It was flrst found experimentally 
|17Llla.ll9| and then investigated theoretically using sim- 
ple models I20I [2l| as well as simulations of the full 
fluid equations |22l |23| | . This fascinating state of spatio- 
temporal chaos is characterized by rotating spirals with 
varying numbers of arms and of different size, which ap- 
pear and disappear irregularly and interact with each 
other and with other defects. SDC arises from the roll 
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state at a threshold that can be as low as e = 0.1 m the 
limit of small Pr. It is driven by large-scale flows that 
are induced by roll curvature and have a strength that is 
proportional to Pr~^ [2^ . 

So far, strongly non-linear NB convection has been 
studied mostly for large Prandtl numbers 0, |3j l25l |. 
but little is known for small ( Pr ~ 1) or very small 
Prandtl numbers (Pr <gC 1). In particular, whether reen- 
trant hexagons exist at large e in the presence of the 
large-scale flows that arise at low Pr, and how NB ef- 
fects impact spiral defect chaos are interesting questions, 
which we address in this paper. 

Here we study NB convection in gases with small 
Prandtl numbers. Specifically, we consider parameters 
corresponding to convection in CO2 and SFg {Pr ~ 0.8) 
and in a H2-Xe mixture {Pr = 0.17). We show that reen- 
trant hexagons are possible in convection in CO2. For 
spiral defect chaos in SFg we find that NB effects pro- 
mote small convection cells ('bubbles'). In H2-Xe, closer 
to threshold, roll-like structures dominate. In both cases 
the NB effects reduce the spiral character of the pattern. 
We quantify the impact of the NB effects on spiral defect 
chaos using geometric diagnostics that we have proposed 
recently |2q|. 

The paper is organized as follows. In Sec^jwc briefly 
review the basic equations, emphasizing how our com- 
putations focus on weakly non-Boussinesq effects, but 
strongly nonlinear convection. In Sec lllll we present the 
results for the linear stability of hexagons and rolls in 
CO2 for a range of parameters accessible experimentally. 
To compare with experiments, we study the influence of 
different lateral boundary conditions on the transition 
from hexagons to rolls in Sec. IIVI In Sec. 0we discuss 
spiral defect chaos in SFg under NB conditions. The sta- 
bility of hexagons and spiral defect chaos in fluids with 
very low Prandlt number {Pr — 0.17) is studied in a mix- 
ture of H2 and Xe in Sec. IVII Finally, conclusions are 
drawn in Sec lVIII 



II. BASIC EQUATIONS 

The basic equations that we use for the description of 
NB convection have been discussed in detail previously 
[16112^ . We therefore give here only a brief summary. We 
consider a horizontal fluid layer of thickness d, density 
p, kinematic viscosity v, heat conductivity A, thermal 
diffusivity k, and specific heat Cp. The system is heated 
from below (at temperature Ti) and cooled from above 
(at temperature T2 < Ti). 

To render the governing equations and boundary con- 
ditions dimensionless we choose the length d, the time 
(P/ko, the velocity Ko/d, the pressure poVoKq/cP, and the 
temperature Tg = voKq/ a^gcP as the respective scales. 
The subscripted quantities refer to the respective values 
at the middle of the fluid layer in the conductive state. 
The non-dimensionalization gives rise to two dimension- 
less quantities: the Prandtl number Pr = i/q/kq, and the 



Rayleigh number R = agATgd^ /voKq. Furthermore, we 
write the equations in terms of the dimensionless momen- 
tum density Vi — pdui/ pqKq instead of the velocities Ui. 
The dimensionless form of the temperature T — T/Tg, 
heat conductivity A = A/Aq, density p ~ p/ Po^ kinematic 
viscosity — vjvQ, and specific heat Cp = CpjcpQ will 
be used in the ensuing equations and the hats dropped 
for clarity. In dimensionless form the equations for the 
momentum, mass conservation and heat are then given, 
respectively, by 
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with the dimensionless boundary conditions 



v{x, y, z, t) = Q{x, y, z, t) = at z = ± 



1 



(4) 



Here <d is the deviation of the temperature field from 
the basic conductive profile. Summation over repeated 
indices is assumed. 

We consider the NB effects to be weak and retain 
in a Taylor expansion of all material properties only 
the leading-order temperature dependence beyond the 
OB approximation. For the density this implies also a 
quadratic term with coefficient 71. It contributes, how- 
ever, only to the buoyancy term in ||2Jl; in all other ex- 
pressions it would constitute only a quadratic correction 
to the leading-order NB effect. Thus, the remaining tem- 
perature dependence of the fluid parameters p, v, A, and 
Cp in (|2I2I3|I is taken to be linear 
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The coefficients 7^ give the difference of the respective 
fluid properties across the layer. They depend therefore 
linearly on the Rayleigh number. 



7.(Ar)=7,= 



7r(l + e), 



(9) 
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where is the value of 7^ at the onset of convection and 
e = {R— Rc{'^f)) / Rc{jf) is the reduced Rayleigh number. 

In analogy to 0, we further omit NB terms that 
contain cubic nonlinearities in the amplitudes Vi or 
Q, as they arise from the expansion of the advection 
terms Vjdj{vi/p) and {vj/p)djQ when the temperature- 
dependence of the density is taken into account. Since 
we will be considering Rayleigh numbers up to twice 
the critical value, which implies enhanced NB effects, 
these approximations may lead to quantitative differ- 
ences compared to the fully NB system, even though the 
temperature-dependence of the material properties them- 
selves is in most situations well described by the linear 
(or quadratic in the case of the density) approximation. 

To quantify the overall strength of the NB-effects we 
use Busse's NB parameter Q, which is given by 

4 

4=0 

where the quantities Vi are linear functions of Pr^^ . The 
NB parameter Q quantifies the breaking of the up-down 
symmetry, which renders at most one of the two types 
of hexagons stable. Gases have a positive value of Q 
and exhibit hexagons with down-flow in the center {g- 
hexagons), whereas liquids have negative Q and show 
hexagons with up-flow (/-hexagons). 

We focus in this paper on the stability properties of 
patterns in the strongly nonlinear regime. They are de- 
termined by a Galerkin method (e.g. [13 )■ We use a 
Fourier expansion on a hexagonal lattice in the lateral 
directions. The Fourier wave vectors q are constructed 
as linear combinations of the hexagon basis vectors bi = 
(7(1, 0) and b2 = 9(1/2, ^3/2) with q — mhi-\-nh2 where 
the integers m and n are in the range |mbi -|-7ib2| < Uqq. 
The largest wavenumber is then Uqq and the number of 
Fourier modes retained is given by 1 -f 6^^^^ j. Typ- 
ically we use Ug = 3. The top and bottom boundary 
conditions are satisfied by using appropriate combina- 
tions of trigonometric and Chandrasekhar functions in z 
(0) Ell)- In most of the computations we use — Q 
modes for each field in Ea. H2l2l3|l . The linear analysis 
yields the critical Rayleigh number Rc as well as the crit- 
ical wavenumber q^. Both depend on the NB-coefficients 
which in turn depend on Rc. Thus, in principle one 
obtains an implicit equation for the 7^^^. The shift in the 
critical Rayleigh number away from the classical value 
Rc = 1708 due to the NB-effects is, however, quite small 
(less than 1 percent) and therefore the resulting change 
in the 7^ is negligible. In this paper we therefore choose 
the 7,f corresponding to Rc = 1708. 

To investigate the nonlinear hexagon solutions, we 
start with the standard weakly nonlinear analysis to de- 
termine the coefficients of the three coupled amplitude 
equations for the modes making up the hexagon pattern. 
To obtain the fully nonlinear solutions requires the solu- 
tion of a set of nonlinear algebraic equations for the ex- 
pansion coefficients with respect to the Galerkin modes. 



To 


AT, Pr 7S 7? 7| 7I 7I Q 


20 


9.43 0.87 0.0486 -0.0669 0.0779 0.0236 -0.0251 1.199 


40 


15.52 0.84 0.0685 -0.0883 0.1132 0.0508 -0.0184 1.724 


60 


23.80 0.82 0.0931 -0.1148 0.1566 0.0919 -0.0074 2.430 



Table I: Values for the critical temperature (in °C), the 
Prandtl number Pr, NB coefficients 7^ , and Busse's parame- 
ter Q for CO2 at the onset of convection for three values of 
the mean temperature (in °C). The values correspond to a 
layer thickness of d = 0.08 cm and a pressure of P = 300 psi. 



This is achieved with a Newton solver for which the 
weakly nonlinear solutions serve as convenient starting 
solutions. In the Galerkin code amplitude instabilities 
are tested by linear perturbations of the expansion coef- 
ficients. In addition, modulational instabilities are con- 
sidered, which involve the introduction of Floquet multi- 
pliers exp(zs • {x,y)) in the Fourier ansatz for the linear 
perturbations of the Galerkin solutions. 

We also study the temporal evolution of the system. 
For that we employ a Fourier spectral code on a rectangu- 
lar grid (idqo-Jdqy), i,j = I...N, with dqy/dq^ = \/3/2 
to accomodate perfect hexagonal patterns. The same ver- 
tical modes are used as in the Galerkin stability code 
|2^ . To solve for the time dependence a fully implicit 
scheme is used for the linear terms, whereas the non- 
linear parts are treated explicitly (second order Adams- 
Bashforth method). The time step is typically taken to 
be ty/500, where ty is the vertical diffusion time. We 
have tested that the stability regimes obtained from the 
Galerkin analysis are consistent with the direct numeri- 
cal simulations. Both codes employed in this paper have 
been kindly provided by W. Pesch . 

III. REENTRANT HEXAGONS IN CO2 

In this paper we investigate specific scenarios for con- 
vection in gases that should be experimentally realizable. 
We focus in this section on CO2 in a conventional range 
for the layer thickness, pressure, and temperature. Table 
n provides the NB coefficients and the Q-value at the on- 
set of convection for a representative range of the mean 
temperature Tq in a layer of thickness d = 0.08 cm [45| . 



A. Amplitude instabilities 

In our analysis we first concentrate on spatially peri- 
odic solutions with their wavenumber being fixed at the 
critical wavenumber and discuss their domains of exis- 
tence and their stability with respect to amplitude insta- 
bilities, which do not change the wavenumber. For three 
different cells with thicknesses d — 0.07, 0.08, and 0.09 
cm, respectively, we consider the range < e < 1 at a 
pressure of P = 300 psi. 
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Figure 1: Stability regions for hexagons and rolls in CO2 with 
respect to amplitude perturbations for three fluid depths: d — 
0.07 cm (dotted lines), d = 0.08 cm (fuU lines), d = 0.09 cm 
(dot-dashed line). Pressure is kept at P = 300 psi. Thick 
lines: stability boundaries for hexagons. Thin lines: stability 
boundaries for rolls. For a given depth, rolls are stable above 
the thin line, and hexagons are unstable in the inner region 
of the thick line. 

The results of the stability analysis for hexagons and 
rolls are shown in Fig^ The hexagons are linearly stable 
for very small e. For a given layer thickness and not too 
high mean temperature Tq the hexagons become unsta- 
ble as the control parameter is increased. The hexagon 
patterns then undergo a second steady bifurcation as the 
control parameter is increased further and become sta- 
ble again. Such restabilized hexagons have been termed 
reentrant hexagons [3 HE US • As the mean tempera- 
ture is increased or the layer thickness is decreased the 
critical heating and with it the NB effects increase. This 
shifts the point of reentrance to lower e and the lower 
stability limit to higher e, decreasing the e-range over 
which the hexagons are unstable, until the two limits 
merge at a temperature Tm- For Tq > the hexagons 
are amplitude-stable over the whole range of e considered 
(0 < e < 1). 



We have also computed the stability of rolls with re- 
spect to amplitude perturbations. The corresponding 
stability limits are indicated in Fig^by thin lines. Rolls 
are stable above these lines. As the NB effects become 
stronger the stabilization of rolls is shifted to larger e. In 
contrast to the hexagons, the rolls do not undergo a sec- 
ond bifurcation within the parameter regime investigated 
and remain amplitude-stable beyond e = 1.0. For strong 
NB effects one has therefore a large range of parameters 
over which the competing rolls and hexagons are both 
linearly amplitude-stable. 

The amplitude-stability limits of the hexagons and 
rolls depend on their wavenumber. This is illustrated 
for the hexagons in Fig. [21 for a mean temperature of 
Tq = 40 °C. The instability region with respect to am- 
plitude perturbations forms a bubble-like closed curve, 
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Figure 2: Stability regions for hexagons in CO2 with respect 
to amplitude (dashed line) and side-band perturbations (solid 
line) for layer depth d = 0.08 cm, pressure P = 300 psi, mean 
temperature To = 40 °C, and Prandtl number Pr = 0.84. The 
NB coeflicients 7f are reported in Tab. Q Hexagons are stable 
with respect to amplitude perturbations outside the dashed- 
line region, and stable with respect to side-band perturbation 
inside the solid-line region. 

inside of which the hexagons are unstable with respect 
to amplitude perturbations. 

It is worth mentioning that the stability limits for 
hexagons in CO2 are quite similar to those of NB convec- 
tion in water, except that in CO2 the NB effects increase 
rather than decrease with increasing mean temperature 

E 



B. Side-Band Instabilities 

In systems with sufficently large aspect ratio side-band 
instabilities can be the most relevant instabilities. Using 
the Galcrkin method, we have studied the stability of the 
hexagons with respect to long- and short-wave side-band 
perturbations for Tq ^ 40°C. The results are shown in 
Fig|21 We find that over the whole range < e < 1 
the only relevant side-band perturbations are long-wave 
and steady, as is the case in the weakly nonlinear regime. 
The same is true also in water with higher Prandtl num- 
ber |l6j| . In this parameter regime the stability region 
consists of two disconnected domains, reflecting the reen- 
trant nature of the hexagons. The stability domain near 
onset is very small and closes up as the amplitude sta- 
bility limit is reached. In the reentrant regime the stable 
domain opens up again in an analogous fashion when the 
amplitude-stability limit is crossed. Note that the sta- 
bility boundaries are leaning toward lower wavenumbers. 
Thus, stable reentrant hexagon patterns are expected to 
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Figure 3: Stability regions for hexagons in CO2 with respect 
to side-band perturbations for layer thickness d — 0.052 cm, 
pressure P = 300 psi, mean temperature To = 27°C, and 
Prandtl Pr = 0.86. The NB coefficients are 7S = 0.2053, 7? = 
-0.2877, 7§ = 0.3267, 7I = 0.1152, and 7I = -0.086. The 
hexagons are stable with respect to side-band perturbations 
in the region inside the solid lines and unstable outside. The 
dashed line corresponds to the neutral curve. 



have wavenumbers below qc- This is related to the fact 
that in the OB case hexagons can be stable for large e, 
but they are side-band stable only for small wavenumbers 

As the mean temperature is increased the bubble of the 
amplitude instability shrinks and eventually disappears, 
as shown in Fig. Ofor a cell with thickness d = 0.052 cm 
and mean temperature Tq = 27 °C. As before, the rele- 
vant side-band instabilities are long-wave all along the 
stability limits and with increasing e the wavenumber 
range over which the hexagons are stable shifts towards 
smaller wavenumbers. For these parameters the region 
of side-band-stable hexagons reaches without interrup- 
tion from the strongly nonlinear regime all the way down 
to threshold. For yet stronger NB effects the range of 
stable wavenumbers widens. 



IV. COMPARISON WITH EXPERIMENTS IN 

CO2 

Bodenschatz et al. carried out a set of experiments 
on convection in CO2 in a cylindrical cell with aspect 
ratio F ~ 172, thickness d = 0.052 cm, and pressure 
P = 300 psi. Under these conditions NB effects are rel- 
evant. In the experiments a weakly hysteretic transition 
from hexagons to rolls was found near e = 0.1. Noting 
that this transition point was below the amplitude insta- 
bility of hexagons to rolls as predicted by weakly nonlin- 
ear theory, the authors interpreted their results in terms 



of the heterogeneous nuclcation of rolls by the sidewalls. 
They found that for small e the concentric rolls induced 
by the sidewall heating remained confined to the inmedi- 
ate vicinity of the sidewalls; however, as e was increased 
the rolls invaded the hexagons and filled the whole cell, 
inducing a transition from hexagons to rolls. 

A comparison of the experimental findings with the 
stability results shown in Fig. O shows that indeed the 
transition cannot be due to an amplitude instability of 
the hexagons. In fact, in this regime the NB effects are so 
strong that, in contrast to the predictions of the weakly 
nonlinear theory, the hexagons do not undergo an am- 
plitude instability at all. To clarify the influence of the 
sidewalls and to assess the significance of the side-band 
instabilities for the transition from hexagons to rolls, we 
perform direct simulations of the Navier-Stokes equations 
II2I2I3I4II for two different sets of boundary conditions 

m 

i) periodic boundary conditions, 

ii) concentric rolls as boundary conditions. In our com- 
putations this type of boundary condition is generated by 
a suitably patterned heating in the interior of the fluid. 
In the experiments concentric rolls near the side walls 
were generated by a side- wall heating }l2i |. 

i) according to Fig. O for e = 0.3 hexagons with 
wavenumbers larger than q > 2.98, which includes the 
critical wavenumber, are unstable with respect to side- 
band instabilities. To test whether these side-band in- 
stabilities trigger a transition to rolls we perform nu- 
merical simulations with periodic boundary conditions 
and hexagons as initial conditions. Fig. ^presents some 
snapshots of the ensuing temporal evoluation in a cell of 
size L — %-2'n/qc — 16.11 using 128 x 128 Fourier modes. 
More precisely, to allow perfect hexagons in a rectangular 
container the container cannot be square. In our simu- 
lations we use Lx = L and Ly = ^/3L/2. The sideband 
instability of the initially almost perfect hexagon pattern 
(cf. Fig. ^a)) induces a shearing of the pattern (cf. 
Fig. ^b)). At the same time a few penta-hepta defects 
arise and some hexagonal convection cells tend to con- 
nect with each other forming short roll-like structures. 
However, as time evolves the system becomes progres- 
sively more ordered again and eventually after losing a 
number of convection cells a defect-free hexagon pattern 
with a smaller wavenumber is established (cf. Fig. QJc) 
at t ~ 60ty). 

Thus, while roll-like features appear for this value of e, 
with periodic boundary conditions no transition to rolls 
occurs and the system relaxes to a new ordered hexagon 
pattern. Only for yet larger values of e the roll-like struc- 
tures that arise in the intermediate, disordered state take 
over and lead to a transition to rolls induced by the side- 
band instabilities. 



ii) clearly, the simulations for periodic boundary con- 
ditions do not match the experimental results described 
above, where a transition from hexagons to rolls occurs 




Figure 4: Numerical simulation corresponding to a layer of 
CO2 with thickness d = 0.052 cm, pressure P = 316 psi, mean 
temperature To = 27.3°C, Prandtl number Pr = 0.87, and 
control parameter e — 0.3 (cf. Fig^J. The size of the inte- 
gration domain is L = 8 ■ 2Ti/qc — 16.11 and the boundary 
conditions are periodic. As initial condition a perfect hexagon 
lattice with wavenumber qc = 3.12 has been used, a) corre- 
sponds to t = 0, b) to f = 18tv, and c) to t — GOtv 



already for e > 0.11. To address this disagreement we 
take into account the fact that in the experiments the 
side walls promote the nucleation of rolls JJl . 

Strictly speaking, the code we are using does not allow 
non-periodic boundary conditions. To mimic the experi- 
mentally used clylindrical cell we employ a step ramp in 
e that reduces e to values well below e = outside a circle 
of radius r = 0.45L with L = 16 • 27r/gc = 32.22 To 
induce concentric rolls near the sidewalls we introduce 
for r > 0.451/ an additional heating in the interior of the 
fluid in the form of concentric rings. Using hexagonal 
initial conditions in the bulk, this leads to an initial state 
as shown in Fig. 

Fig. [SJi,b shows two snapshots at t — ty and t = 158ty 
demonstrating how the rolls induced by the side walls in- 
vade the carefully prepared hexagonal state in the bulk 
already for e = 0.2. This is well below the e- value for 
which with periodic boundary conditions the hexagons 
persisted even through the side-band instability. The fi- 
nal steady state consists of concentric rolls as observed 
in the experiments (cf. Fig. 5 in [l^'). For lower e, how- 
ever, the experimentally observed final state consists of 
hexagons in the bulk of the system surrounded by wall- 
induced concentric rolls (cf. Fig. 4 in 12]). We find 
this also in our numerical simulations, as shown in Fig. 
Efc-d). There the forcing of rolls is identical to that in 
Fig. IHl^a-b) but e — 0.05. Starting with random initial 
conditions, the forcing gives rise to a ring contained in 
the square integration domain. At the beginnig of the 
simulation (Fig. EIc)) the rolls created by the forcing in- 
vade the interior of this small system. However, as time 
progresses the rolls pull out of the central region of the 
cell, and the final steady state (for t > 165 ti,) consists 
of stable hexagons surrounded by a couple of concentric 
rolls in addition to those induced by the forcing. 



Thus, our simulations suggest that the experimentally 
observed transition from hexagons to rolls is neither due 
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Figure 5; Numerical simulation in a cell of CO2 with thick- 
ness d — 0.052 cm, pressure P — 300 psi, mean temperature 
To = 27°C, and Prandtl number Pr = 0.86. The size of the 
integration domain is L = 16 ■ 2 vr/gc ~ 32.22, and the bound- 
ary conditions are rings of concentric rolls generated with and 
external forcing. 

For e = 0.2 (upper snapshots) hexagon inititial conditions 
have been used, with t — (left snapshot) and t — 158 1„ 
(right). The lower snapshots correspond to a simulation with 
e — 0.05 and random initial conditions, at t = 40 tv (left) and 
t = 158 1„ (right). 



to amplitude instabilities nor to side-band instabilities. 
Rather, there is a large range of parameters in which 
hexagons and rolls are both linearly stable and the final 
state is selected by one type of pattern invading the other. 
The transition to rolls at these low values of e is made 
possible by the boundaries, which provide a seed for the 
nucleation of rolls. We expect that by applying a forcing 
that is confined to the region near the boundaries and 
that replaces the rolls by hexagons the transition to rolls 
could be shifted to substantially larger values of e. Such 
a forcing could be achieved by a patterned heating of 
the interior of the fluid jsj or possibly by a suitable 
geometric patterning of the bottom plate |3^ . 



V. NB SPIRAL DEFECT CHAOS IN SFe 

A fascinating state observed in convection at low 
Prandtl numbers is spiral defect chaos. It is character- 
ized by self-sustained chaotic dynamics of rotating spi- 
rals, as well as dislocations and disclinations, and arises 
in fluids with Pr<l[ll[ll|2l|2l|2lina parame- 
ter range where straight rolls are linearly stable. Spiral 
defect chaos has so far predominantly been investigated 
under OB conditions, in which up-flows and down-flows 
are equivalent. 

As mentioned before, NB effects break the up-down 
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symmetry and different flow structures may be predom- 
inantly associated with up-flow and down-flow, respec- 
tively. Moreover, in the absence of the OB symmetry 
a resonant triad interaction is allowed. If it is strong 
enough it leads to the formation of hexagons. For weaker 
interaction one may still expect an enhancement of cel- 
lular rather than roll-like structures. 

To investigate the impact of NB effects on spiral de- 
fect chaos we consider convection in a layer of SFg. This 
gas has been used previously in experimental convection 
studies under usual laboratory conditions [33|. and near 
the thermodynamical critical point 14. ilfl l84l|. In Fig. 
El^a) we present the stability diagram for hexagons and 
rolls with respect to amplitude perturbations in a layer 
of SFg of thickness d = 0.0542 cm, pressure P — 140 psi, 
and a range of temperatures that is experimentally ac- 
cessible. Hexagons are amplitude-stable to the right of 
the solid line and rolls above the dashed line. As in the 
case of CO2 the NB effects increase with increasing mean 
temperature Tq and above a certain value of Tq hexagons 
are linearly amplitude-stable over the whole range of e 
investigated. Here we focus on relatively strong NB ef- 
fects. We therefore show in Fig|Hb the stability limits 
with respect to side-band perturbations for a relatively 
large mean temperature, Tq — 80°C. As in the case of 
CO2, the wavenumber range over which the hexagons are 
stable is leaning towards smaller wavenumbers. Overall, 
amplitude and side-band stability limits are qualitatively 
similar to those of convection in CO2 (cf. Fig. P). 



Fig[7| shows two snapshots obtained by direct numer- 
ical simulations of the Navier-Stokes equations corre- 
sponding to convection in SFg for Tq = 80 °C and e — lA 
in a convective cell of thickness d = 0.0542 cm and hor- 
izontal size L = 16 • 2tt /qc ~ 32.22. Periodic boundary 
conditions are used with 128 x 128 Fourier modes and 6 
vertical modes. Both states are obtained after an inte- 
gration time of 160 ti,, starting from random initial con- 
ditions. While in FiglTja all NB are retained, in Fig[7^ 
the same values are used for Pr and e, but all NB pa- 
rameters 7f are set to 0, i.e. the system is treated as if 
it was Boussinesq. 



The snapshots depicted in Fig[7| show, as expected, 
that due to the NB effects down- flow convection cells, 
which are white in FigEI outnumber cells with up-flow 
(black) . Moreover, in this regime the NB effects enhance 
the overall cellular rather than roll-like character of SDC. 
This manifests itself in the appearance of numerous small 
down-flow convection cells (white 'bubbles') and in the 
appearance of quite noticeable bulges on the NB convec- 
tion rolls. To quantify these and other differences we 
analyse a long sequence of snapshots with a recently in- 
troduced geometric approach '26^. It is based on the 
contour lines corresponding to the intensity half-way be- 
tween the minimal and maximal intensity of all snapshots 




Mean Temperature 




2.4 2.6 2.8 3 3.2 

Wavenumber q 



Figure 6: Stability regions of hexagons and rolls in a layer 
of SFe of thickness d = 0.0542 cm, pressure P — 140 psi, and 
Prandtl number Pr = 0.8. 

a) Stability regions with respect to amplitude perturbations. 
Continous line: stability boundary for hexagons. Dashed line: 
stability boundary for rolls. Stability limits obtained for the 
critical wavenumber Qc. 

b) Stability regions with respect to side-band perturbations 
for the above layer with a mean temperature To — 80 °C The 
corresponding NB coefficients are 70 — 0.1714, 71 — —0.2118, 
72 = 0.2836, 7i = 0.1905, and 74 = 0.0624 corresponding to 
Q = 4.2.. The dashed line corresponds to the neutral curve. 

in a given run. The contour lines corresponding to the 
temperature field of snapshots Fig[7Ja,b) are shown in 
Fig[7Ic,d). In the following we present various statistics 
of these contour lines. 



The most striking difference between the OB and the 
NB case is the asymmetry that is induced by the NB 
effects between 'black' and 'white' components, i.e. be- 
tween closed contours that enclose up- and down-flow re- 
giones, respectively. To quantify this asymmetry we mea- 
sure the number of white and black components. These 
two topological measures correspond, respectively, to the 
Betti numbers of order 1 and 2 of the pattern defined by 
the white components [s^. Fig|Sl shows these two quan- 
tities as a function of time in the OB and the NB case. 
As expected, in the OB case the number of black and 
white components is essentially the same at all times, 
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OB NB 




Figure 7: Direct numerical simulation of equations I2I4II for 
SFe in a cell of thickness d = 0.0542 cm, pressure P = 140 psi, 
mean temperature To = 80°C, Prandtl number Pr — 0.8, 
and control parameter e = 1.4. The cell has size L — 16 ■ 
2n/qc = 32.22 with periodic boundary conditions. Starting 
from random initial conditions both snapshots are taken at 
t = 160 tv. Left panels for OB conditions (7i=0, i = 0..4, left), 
right panels for NB conditions appropriate for To = 80° C (cf. 
Fi^Sja). Bottom panels give the corresponding contour lines 
used for the pattern diagnostics. 




Time (t ) 

Figure 8: Number of black and white closed contours as a 
function of time for NB and OB conditions for the simulations 
corresponding to the snapshots shown in Fig. 




0.2 0.4 0.6 0.8 1 1.2 1.4 
Red. Rayleigh Number 8 

Figure 9: Total number of closed contours as a function of 
the control parameter e for NB (circles) and OB (squares) 
conditions. 

interaction that is made possible by the breaking of the 
OB symmetry and which tends to enhance cellular rather 
than filamentary roll-like structures. 



whereas in the NB case the white components signifi- 
cantly outnumber the black ones. The ratio of white to 
black components is therefore a sensitive indicator for 
the significance of NB effects. FiglHlalso illustrates how 
much the number of components fluctuates during these 
runs. Recently, the two Betti numbers have also been 
measured based on patterns obtained in experiments on 
SDC in convection in CO2. Scanning e in very small steps 
the authors report steps in the Betti numbers indicative 
of transitions between different chaotic, disordered states 




FiglHl shows that the total number of components 
(closed contours) is considerably larger in the NB case 
than in the OB case. This is presented in more quan- 
titative detail in Fig|21 which gives the mean value of 
the total number of components, i.e. the sum of black 
and white components, as a function of e for OB as well 
as NB conditions. In the NB case the total number of 
components is up to four times larger than in the OB 
case. We attribute this difference to the resonant triad 



To characterize the components better and to distin- 
guish cellular and roll-like structures we introduce the 
'compactness' C of components p^ . 

C = 47r^- (11) 

Here A is the area inside a closed contour and P its 
perimeter. With the normalization used in (|ll|l compact, 
cellular structures are charecterized by C < 1, whereas 
filamentary, roll- like structures have C <^1. 



Fig. 1 101 shows the mean number of closed contours per 
snapshot for a given compactness C for the NB and the 
OB simulation at e = 1.4 over the duration t = 360 1^,. 
As expected, in the NB case the number of white com- 
ponents is much larger than that of black components, 
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Figure 10: Mean number of closed contours per snapshot 
with a given compactness C for convection in SFe at e = 1.4 
under NB (circles) and under OB (squares) conditions (cf. 

Fig.m. 
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Figure 11: Distribution function for the compactness C for 
NB (circles) and OB (squares) conditions (cf. Fig. |7J. 



whereas in the OB case both are about the same. The 
total number of components is noticeably larger in the 
NB case, which also shows an increase in the number of 
white, filamentary contours with small compactness. 



The relative distribution among components with dif- 
ferent compactness is more clearly visible in the relative 
frequency of contours as a function of the compactness, 
which is shown in Fig^] More precisely, its data re- 
sult from running averages over adjacent 10 points of the 
corresponding data shown in Fig llOl which are then nor- 
malized. The normalized data show that the increase in 
the relative frequency of white filamentary components 
[C ^ 1) is essentially the same in the NB case and in 
the OB case. However, only very few black filamentary 
components arise in the NB case. 

A feature of Fig^Jthat is surprising at first sight is the 
essentially equal height of the NB peak and the OB peak 
for components with C ^ 1. Visually, the NB snapshot 
exhibits many more small compact 'bubbles' than the 



OB run. An explanation for this observation can be ob- 
tained by correlating the compactness of the closed con- 
tours with their length V . The joint distribution function 
for these two quantities is shown in Fig ll2l using logarith- 
mic scales for V and C. Focussing on the compact objects 
with C < 1 one recognizes in the OB case a second peak at 
somewhat larger contourlength. We associate this peak 
with the appearance of target-like structures, i.e. with a 
second contour line encircling a smaller compact, almost 
circular contourline. In the NB case this second peak is 
barely visible. Instead, the shoulder of the main peak 
is extended significantly towards smaller contourlength 
V. It signifies the appearance of compact objects that 
are smaller than a typical wavelength, which we asso- 
ciate with the small 'bubbles' that are easily recognized 
in the snapshot of the NB case Fig. ^}p. Thus, the com- 
parable relative frequency of small components shown in 
FigEJhas a different origin in the OB and the NB case. 
Whereas in the NB case it is mostly due to small bub- 
bles, it seems to originate from target-like structures in 
the OB case. 

In the logarithmic scaling used in FigEl a straight 
ridge arises in the distribution function for large con- 
tourlengths. It is characteristic for filamentary structures 
with a typical width, which corresponds here to half a 
wavelength A of the convection rolls. 



C 



(12) 



In the OB case one can discern deviations from this scal- 
ing. They are confined to larger rather than smaller com- 
pactness values for a given contourlength, indicating that 
the long rolls can be wider but not narrower than a cer- 
tain thickness. In the NB case these deviations are much 
smaller. 

To identify spiral components in the pattern directly 
we also measure the winding number of the components 
[2^ . It is defined via the angle 9 by which the (spiral) 
arm of a pattern component is rotated from its tip to 
its end at the vertex at which it merges with the rest 
of the component. In cases in which a component has 
no vertices we split it into two arms at the location of 
minimal curvature 26] . The winding number is then de- 
fined as \W\ = 9/2TT. To assess the impact of the NB 
effects on the spiral character of the pattern we measure 
the number of spirals in each snapshot and show the re- 
sulting histogram over the whole run in Fig ll3l We use 
three different thresholds Wmin for the identification of 
spirals, Wmin = 1, >V,„i„ = 1/2, and Wmin = 1/4. As 
Figll3lshows. the number of small spirals with |>V| > 1/4 
is quite similar in the OB and the NB case. However, 
larger spirals with |W| > 1/2 or even |>V| > 1 are much 
more rare in the NB case; in fact, for the system size 
L = 16 ■ 27r/qc = 32.22 that we have used in these simu- 
lations there was at most one spiral with \ W\ > 1 at any 
given time. 

The reduced spiral character of NB spiral defect chaos 
is also quite apparent in Fig ll4l which shows the cor- 
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Figure 12: Distribution function for contour length and com- 
pactness for OB (a) and NB conditions (b) Pr — 0.8 and 
e = 1.4. 
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Figure 14: Correlation between the arclength and the wind- 
ing number of spiral arms under OB (black squares) and NB 
conditions (red circles). 



relation between the winding number and the arclength 
of the spiral arm. More precisely, each dot marks the 
occurrence of one spiral arm in a snapshot. In the OB 
case one can see quite clearly a maximal winding num- 
ber for any given arclength, which is consistent with an 
Archimedean shape of the spiral 26J 47J. In the NB case, 
however, only components with very small contourlength 
reach the Archimedean limit and most components have 
winding numbers that are much smaller, i.e. the compo- 
nents are quite straight. 



VI. HEXAGONS AND SPIRAL DEFECT 
CHAOS AT VERY LOW PRANDTL NUMBERS. 
Ha-XE MIXTURES 
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Figure 13: Distribution function for the number of spirals 
in spiral defect chaos under OB and NB conditions for three 
values of the threshold, Wmin = 1/4 (circles), Wmin = 1/2 
(squares), and Wmin = 1 (diamonds). 



As mentioned above, the restabilization of NB 
hexagons at larger Rayleigh number is related to the ex- 
istence of stable OB hexagons at large Rayleigh numbers. 
The wavenumber range over which the OB hexagons are 
stable shrinks with decreasing Prandtl number and for 
Pr < 1.2 the OB hexagons are side-band unstable at all 
wavenumbers j3Qi] . However, as seen in the case of CO2 
and SFg, NB hexagons can be side-band stable at large 
e even below Pr = 1.2 due to the additional stabilizing 
effect of the resonant triad interaction. It is therefore of 
interest to investigate whether the NB effects can be suf- 
ficient to stabilize strongly nonlinear hexagons even for 
Prandtl numbers significantly below Pr = 1. 

Prandtl numbers well below Pr — 1 can be reached by 
using a mixture of a heavy and a light gas. An experi- 
mentally investigated case is a mixture of H2 and Xe [s^l ■ 
With a mole fraction of x = 0.4, one can reach Prandtl 
numbers as small as Pr — 0.17. The Lewis number of 
such a mixture is close to one 0, . Therefore such 
mixtures are expected to behave essentially like a pure 
fluid with the same Prandtl number 

We investigate the stability of hexagon and roll con- 
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vection in a H2-Xe mixture with mole fraction % = 0.4 at 
a pressure of 300 psi and a layer thickness of d — 0.1cm. 
With respect to amplitude instabilities the stability di- 
agram is very similar to that of convection in CO2 and 
SFg with hexagons becoming reentrant at e-values as low 
as e = 0.14 for Tq = 20 °C. 

Focusing on strong NB effects we perform a detailed 
stability analysis with respect to sideband perturbations 
at a mean temperature of Tq — 80°C using the same 
layer thickness oi d = 0.1cm. For these low Prandtl 
numbers the numerical resolution has to be increased to 
obtain sufficiently well resolved flows. While for the sta- 
bility anlyses of hexagons in CO2 and SFg it is sufficient 
to use Uz = 6 and = 3 in the Galerkin expansion, 
for the H2-Xe mixture with Pr ~ 0.17 at least Ug = 5 
and 71,2 = 6 are required. Fig. depicts the result- 
ing stability diagram. It shows that the region of side- 
band stable hexagons is not contiguous but consists of 
the usual region immediately above threshold and an ad- 
ditional, disconnected region at larger Rayleigh numbers. 
We could not follow the stability limits to smaller values 
of q than shown in Fig llSI due to numerical convergence 
problems. Presumably, these arise due to bifurcations in- 
volving additional, resonant wavevectors jS^ , somewhat 
similar to the resonances studied in Taylor vortex flow 
|40l l4ll |. The fact that the region of stability is dis- 
connected is remarkable since the region of amplitude- 
stability (not shown) is actually contiguous. This is in 
contrast to the behavior found in CO2 and SFg where 
the two side-band stable regions become connected when 
the bubble-like region of amplitude instability disappears 
(cf. Fig l2l3|l . The comparison of the stability limits with 
those of CO2 and of SFe (Figgis) shows further that the 
maximal wavenumber q at which the hexagons are stable 
with respect to side-band perturbations decreases with 
decreasing Prandtl number and as a result the over-all 
stability region shrinks as well. 



The stability analysis of the H2-Xe mixture also re- 
veals an oscillatory instability of the hexagon patterns at 
e ~ 1. Within the e-range investigated, no such oscil- 
latory instability was found at the larger Prandtl num- 
bers relevant for CO2 and SFg. Unfortunately, it turns 
out that before the onset of the oscillatory instability the 
hexagons already become unstable to a side-band insta- 
bility at half the hexagon wavelength, which will always 
preempt the oscillatory instability. For the rolls we also 
find an oscillatory instability. It is presumably related to 
the well-known oscillatory instability of Boussinesq rolls 

m 



For not too small e generic initial conditions will not 
lead to hexagonal patterns but rather to spiral defect 
chaos. A typical snapshot of a pattern in this state is 
shown in Fig^|for e = 0.3. Compared to the patterns 




Wavenumber q 

Figure 15: Stability limits with respect to sideband pertur- 
bations in a mixture of H2 — Xe with molar fraction x = 0.4, 
thickness d — 0.1cm, pressure P = 300 psi, mean tempera- 
ture To = 80 °C, and Prandtl number Pr = 0.17. The NB 
coefficients are 70 = 0.5535, 71 = -0.6421, 72 = 0.9224, 
73 = 0.3647, 74 = -0.0712 resulting in Q = 13.8. Hexagons 
are stable in the regions enclosed by the solid lines and un- 
stable outside. The dashed line represents the neutral curve. 
The results for H2-Xe are obtained with Uz — G and Uq — 5 . 




Figure 16: Typical pattern for convection in H2-Xe at e = 0.3. 



obtained in NB convection in SFg (cf. FigEb) the pat- 
terns in H2-XC arc less cellular and do not show a large 
number of small bubbles. To quantify these and other 
characteristics of the patterns we again apply the geo- 
metric diagnostics introduced earlier l2d| . 

In Fig ll7l we show the normalized distribution func- 
tions for the compactness of white and black components 
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Figure 17; Distribution function of the compactness of closed 
contours for H2-Xe (Pr = 0.17, e = 0.3) . 
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Figure 18: Joint distribution function for the contour length 
and the compactness of closed contours for H2-Xe {Pr = 0.17, 
e = 0.3). . 



(cf. Fig lll|l . Since there are only very few black compo- 
nents their distribution function exhibits large statistical 
fluctuations. Of particular interest is the distribution 
function for the white components. It confirms the vi- 
sual impression that the number of compact components 
is significantly reduced compared to the case of SFg; in 
fact, while in SFg the maximum of the distribution func- 
tion is close to C = 1 , in H2-Xe the absolute maximum is 
at C 0.1, which corresponds to filamentary structures. 
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Figure 19: Distribution function for the absolute value of the 
winding number for NB convection in H2-Xe in comparison 
with convection in SFe in the NB and the OB case. . 



Pi? = 0.17) as well as for Boussinesq and non-Boussinesq 
convection in SFg (at e = 1.4 and Pr = 0.8). As had 
been noted previously in the Boussinesq case 0| the 
distribution function is roughly consistent with exponen- 
tial behavior. In the NB case the exponential decays 
substantially faster than in the Boussinesq case and spi- 
rals with winding numbers above = 1 are rare. In 
the Boussinesq case we had found that the decay rate 
depends mostly on the Prandtl number, but only very 
little on e |2g|. Unfortunately, we do not have enough 
non-Boussinesq data to investigate such trends in the e- 
and Pr-dependence. However, it is worth noting that in 
the two non-Boussinesq cases shown in Fig^l the decay 
rates are essentially the same despite their substantial 
difference in Prandtl number and both decays are much 
faster than that in the Boussinesq case. Thus, possibly 
the impact of NB effects dominates the dependence on 
the Prandtl number. 



In NB convection in H2-Xe, the distribution function 
for the number of spirals is qualitatively similar to that 
in SFg, which is shown in Fig ll3l above. In particu- 
lar, there are almost no spirals with a winding number 
\W\ > 1. Similarly, the correlation between the arclcngth 
and the winding number of spirals reveals that in NB 
convection in H2-Xe there is no significant trend towards 
Archimedean spirals. 



The lack of small bubbles is demonstrated in more 
detail in the joint distribution function for the con- 
tourlength and the compactness, which is shown in 
FigEl Note that the view is rotated compared to FiglT^ 
The distribution function is lacking the broad shoulder 
seen in NB convection in SFg (see Figll2b). Instead, the 
decay of the distribution function towards long filamen- 
tary contours is quite slow. 



To assess the spiral character of NB spiral defect chaos 
at these low Prandtl numbers we show in Fig^|the dis- 
tribution function for the absolute value |>V| of the wind- 
ing number for NB convection in Il2-Xe (at e = 0.3 and 



VII. CONCLUSION 

In this paper we have studied non-Boussinesq convec- 
tion in gases (CO2, SFg, H2-Xe) with Prandtl numbers 
ranging from Pr ~ 1 down to Pr — 0.17 in experimen- 
tally relevant parameter regimes. We have complemented 
a Galerkin stability analysis of hexagon patterns with di- 
rect numerical simulations of the fluid equations to study 
transitions between different hexagon patterns and to 
quantify the impact of non-Boussinesq effects on spiral 
defect chaos. 

We find that the reentrance of hexagons that we have 
identified previously in non-Boussinesq convection in wa- 
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ter [ig also occurs at low Prandtl numbers. As was 
the case at large Prandtl numbers, compressibility is 
not neccessary for reentrance. Since, in addition, the 
range of wavenumbers for which the reentrant hexagons 
are stable differs significantly from that of the reentrant 
hexagons observed in experiments on convection in SFg 
near the thermodynamic critical point the mecha- 
nisms underlying the two types of restabilization of the 
hexagons is most likely different. Reflecting the fact that 
in gases the non-Boussinesq effects increase with increas- 
ing temperature the reentrance is shifted to lower values 
of the Rayleigh number when the mean temperature is 
increased, opposite to the behavior in water [l^ . In con- 
vection in water the reentrant hexagons are stable only 
for wavenumbers below the critical wavenumber. This 
trend becomes more pronounced with decreasing Prandtl 
number. In fact, for the gas mixture with Pr = 0.17 the 
wavenumber at the stability limit of the hexagons de- 
creases so rapidly with increasing Rayleigh number that 
the range in Rayleigh number over which the hexagons 
are stable becomes quite small. 

The comparison of our stability results with experi- 
ments on the transition between hexagons and rolls in 
CO2 [13 shows that this transition is not due to an am- 
plitude or a side-band instability. As a matter of fact, for 
the parameters of the experimental system the hexagons 
do not undergo any linear amplitude instability to rolls, 
contrary to the prediction of the weakly nonlinear the- 
ory |l2l |. We have performed detailed numerical simula- 
tions with various lateral boundary conditions and con- 
firm that the transition is the result of the heterogeneous 
nucleation of rolls at the side walls of the container, which 
then invade the whole system if the Rayleigh number is 
sufficiently high. Our simulations suggest that hexagons 
could be stabilized well beyond the experimentally ob- 
served transition point if the influence of the lateral walls 
can be reduced by applying a spatially patterned forcing 
that drives hexagons at the wall. Such a forcing can be 
achieved by localized heating 'si'l or by geometric pat- 
terning of the bottom plate 32] . Of course, the wavenum- 
ber of the forced hexagons would have to be adjusted to 
lie in the stable range. 

We have also investigated the stability of hexagons in 
H2-Xe mixtures with very small Prandlt number (Pr = 
0.17). There also stable reentrant hexagons are possible, 
but they are restricted to a small range in wavenum- 
ber [q < Qc) and Rayleigh number. Since for such small 
Prandtl numbers Boussinesq hexagons are always side- 
band unstable ^3 j the stability of the hexagons is a re- 
sult of the non-Boussinesq effects. 

A fascinating state that is characteristic for convec- 
tion at low Prandtl numbers is spiral defect chaos [T^ . 
We have studied the influence of NB effects on spiral de- 
fect chaos in a set-up corresponding to convection in SFg 
with a Prandtl number of Pr — 0.8 and in Il2-Xe with 
a Prandtl number of Pr = 0.17. To quantify the dif- 
ferences between Boussinesq and non-Boussinesq spiral 
defect chaos we have employed a recently suggested set 



of geometric diagnostics of the patterns [26j . 

As expected, in SFg and in Il2-Xe the non-Boussinesq 
effects break the equivalence of up- and down-flows and, 
consequently, the mean number of pattern components 
(closed contours) corresponding to down-flow differs sig- 
nificantly from that of the up-flow components. More 
interesting is our finding that in SFg the total number 
of components (Betti number [35l | ) is more than twice as 
large in the non-Boussinesq case than in the Boussinesq 
case for otherwise equal parameters. We attribute this 
enhancement of cellular rather than roll-like structures 
to the resonant triad interaction that is introduced by 
the non-Boussinesq effects. 

Another striking difference between Boussinesq and 
non-Boussinesq spiral defect chaos in SFg becomes ap- 
parent in the joint distribution function for the con- 
tourlength of the components and their compactness. 
While the Boussinesq case exhibits a strong signature of 
target-like structures, the non-Boussinesq case displays 
instead a marked increase in the number of small com- 
pact components ('bubbles'). This trend towards a more 
cellular structure is also apparent in the analysis of spi- 
rals. The number of spirals with a winding number above 
1/2 is much smaller in the non-Boussinesq than in the 
Boussinesq case. In the correlations between the wind- 
ing number and the arclength of the components this is 
reflected by the missing of large Archimedean spirals in 
the non-Boussinesq case. 

Interestingly, we do not flnd many small bubbles in our 
simulations of H2-Xe at Pr — 0.17. In fact, the maximum 
of the distribution function for the compactness of pat- 
tern components is shifted away from C<ltoC~0.1, 
i.e. to fllamentary structures. We do not know whether 
this implies that the tendency towards small bubbles is 
maximal at moderately small Prandtl numbers, i.e. for 
Pr 1, or whether it is due to the fact that for com- 
putational reasons the simulations in H2-Xe were per- 
formed at a lower Rayleigh number (e = 0.3) than in SFg 
{e = lA). 

Due to computational limitations we have not been 
able to investigate the transition from hexagons to spi- 
ral defect chaos at very small Prandtl numbers. Since 
for very small Prandtl numbers there are essentially no 
large spirals, the transition from hexagons to spiral de- 
fect chaos may follow a different path than that observed 
experimentally in CO2 , where large spirals were reported 
to arise in the transition from hexagons to rolls [l3 ■ We 
surmise that the geometric diagnostics that we have em- 
ployed to characterize the well developed spiral defect 
chaos would provide also valuable insight into this tran- 
sition. Since the diagnostics require a substantial amount 
of data the transition would be best investigated experi- 
mentally. 

We deeply appreciate the support by W. Pesch and 
his students, who have developed the codes we have used 
in this study We thank G. Ahlers for providing 

us with the code to determine the NB coefficient. We 
have benefitted from discussions with G. Ahlers, G. Gu- 



14 



naratne, K. Krishan, K. Mischaikow, W. Pesch, and M. FG02-92ER14303) and NSF (DMS-9804673) is gratefully 
Schatz. Support from the Department of Energy (DE- acknowledged. 



[1] E. Bodenschatz, W. Pesch, and G. Ahlers, Ann. Rev. 

Fluid Mech. 32, 709 (2000). 
[2] F. H. Busse, J. Fluid Mech. 30, 625 (1967). 
[3] E. Palm, J. Fluid Mech. 8, 183 (1960). 
[4] L. A. Segel and J. T. Stuart, J. Fluid Mech. 13, 289 

(1962). 

[5] L. A. Segel, J. Fluid Mech. 21, 359 (1965). 
[6] E. Palm, T. EUingsen, and B. Gjevik, J. Fluid Mech. 30, 
651 (1967). 

[7] S. H. Davis and L. A. Segel, Phys. Fluids 11, 470 (1968). 
[8] E. F. Somerscales and T. S. Dougherty, J. Fluid Mech. 

42, 755 (1970). 
[9] M. Dubois, P. Berge, and J. Wesfreid, J. Physique 39, 

1253 (1978). 

[10] F. M. Richter, J. Fluid Mech. 89, 553 (1978). 
[11] R. W. Walden and G. Ahlers, J. Fluid Mech. 109, 89 
(1981). 

[12] E. Bodenschatz, J. R. deBruyn, G. Ahlers, and D. S. 

CanneU, Phys. Rev. Lett. 67, 3078 (1991). 
[13] E. Pampaloni, C. Perez-Garia, L. Albavetti, and S. Cilib- 

erto, J. Fluid Mech. 234, 393 (1992). 
[14] M. Assenheimer and V. Steinberg, Phys. Rev. Lett. 76, 

756 (1996). 

[15] A. Roy and V. Steinberg, Phys. Rev. Lett. 88, 244503 
(2002). 

[16] S. Madruga, H. Riecke, and W. Pesch, J. Fluid Mech. (to 
appear) . 

[17] S. W. Morris, E. Bodenschatz, D. S. CanneU, and 

G. Ahlers, Phys. Rev. Lett. 71, 2026 (1993). 
[18] S. W. Morris, E. Bodenschatz, D. S. CanneU, and 

G. Ahlers, Physica D 97, 164 (1996). 
[19] J. Liu and G. Ahlers, Phys. Rev. Lett. 77, 3126 (1996). 
[20] H. W. Xi, J. D. Gunton, and J. Vinals, Phys. Rev. Lett. 

71, 2030 (1993). 
[21] M. C. Cross and Y. Tu, Phys. Rev. Lett. 75, 834 (1995). 
[22] W. Decker, W. Pesch, and A. Weber, Phys. Rev. Lett. 

73, 648 (1994). 

[23] K.-H. Chiam, M. R. Paul, M. C. Cross, and H. S. Green- 
side, Phys. Rev. E 67, 056206 (2003). 

[24] Y.-N. Young, H. Riecke, and W. Pesch, New J. Phys. 5, 
135 (2003). 

[25] S. Madruga and C. Perez-Garcia, Int. J. Bifurcation 

Chaos 14, 107 (2004). 
[26] H. Riecke and S. Madruga, Chaos (to appear). 



[27] F. H. Busse and R. M. Clever, J. Fluid Mech. 91, 319 
(1979). 

[28] S. Chandrasekhar, Hydrodynamic and Hydromagnetic 

Stability (Clarendon, Oxford, 1961). 
[29] W. Pesch, Chaos 6, 348 (1996). 

[30] R. M. Clever and F. H. Busse, Phys. Rev. E 53, R2037 
(1996). 

[31] D. Semwogerere and M. F. Schatz, Phys. Rev. Lett. 88, 

054501 (2002). 
[32] E. Bodenschatz (priv. comm.). 

[33] E. Bodenschatz, D. CanneU, J. R. DeBruyn, R. Ecke, 
Y. Hu, K. Lerman, and G. Ahlers, Physica D 61, 77 
(1992). 

[34] M. Assenheimer and V. Steinberg, Nature 367, 345 
(1994). 

[35] M. Gameiro, K. Mischaikow, and W. Kalies, Phys. Rev. 
E 70, 035203(R) (2004). 

[36] K. Krishan, M. Gameiro, K. Mischaikow, and M. F. 
Schatz (unpublished). 

[37] J. Liu and G. Ahlers, Phys. Rev. E 55, 6950 (1997). 

[38] G. Ahlers, in Dynamics of spatio-temporal cellular struc- 
tures - Henri Benard centenary review, edited by I. Mu- 
tabazi, J. E. Wesfreid, and E. Guyon (Springer, 2005), 
Springer Tracts in Modern Physics. 

[39] V. V. Moroz, Ph.D. thesis. Northwestern University 
(2004). 

[40] H. Riecke and H.-G. Paap, Phys. Rev. A 33, 547 (1986). 
[41] H.-G. Paap and H. Riecke, Phys. Rev. A 41, 1943 (1990). 
[42] V. Croquette and H. Williams, Phys. Rev. A 39, R2765 
(1989). 

[43] R. Ecke and Y. Hu, Physica A 239, 174 (1997). 
[44] B. Plapp, Ph.D. thesis, CorneU University (1997). 
[45] These values were obtained with a code kindly provided 
by G. Ahlers. 

[46] In the experimental setup the top temperature is held 
constant at T = 12.84°C, and therefore the mean tem- 
perature changes as e is increased. In our computations, 
however, we keep To fixed. Since the transition occurs 
quite close to threshold this is a good aproximation to 
the experimental procedure. 

[47] Note that detailed analyses of large spirals show devia- 
tions from the Archimedean shape due to the dislocations 
that accompany finite spirals |44| 



